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ABSTRACT 

Context, SMM J04542 - 0301 is an extended (~ 1') submm source located near the core of the cluster MS0451.6 - 0305. It has been 
suggested that part of its emission arises from the interaction between a LBG and two EROs at z ~ 2.9 that are multiply-imaged in the 
optical/NIR observations. However, the dramatic resolution difference between the sub-mm map and the optical/NIR images make it 
difficult to confirm this hypothesis. 

Aims. In a previous paper, we reported the detection of 1.4 GHz continuum radio emission coincident with this sub-mm source using 
VLA archival data. To fully understand the relation between this radio emission, the sub-mm emission, and the optical/IR multiply- 
imaged sources, we have re-observed the cluster with the VLA at higher resolution. 

Methods. The previous archival data has been re-reduced and combined with the new observations to produced a deep (~ 10 //Jy 
beam -1 ), high resolution (~ 2") map centred on the cluster core. The strong lensing effect in the radio data has been quantified by 
constructing a new lens model of the cluster. 

Results. From the high resolution map we have robustly identified six radio sources located within SMM J04542 - 0301. The 
brightest and most extended of these sources (RJ) is located in the middle of the sub-mm emission, and has no obvious counterpart in 
the optical/NIR. Three other detections (El, E2 and E3) seem to be associated with the images of one of the EROs (B), although the 
NIR and radio emission appear to originate at slightly different positions in the source plane. The last two detections (CR1 and CR2), 
for which no optical/NIR counterpart have been found, seem to constitute two relatively compact emitting regions embedded in a ~ 5" 
extended radio source located at the position of the sub-mm peak. The presence of this extended component (which contributes 38% 
of the total radio flux in this region) can only be explained if it is being produced by a lensed region of dust obscured star formation 
in the center of the merger. A comparison between the radio and sub-mm data at the same resolution suggests that El, E2, E3, CR1 
and CR2 are associated with the sub-mm emission. 

Conclusions. The radio observations presented in this paper provide strong observational evidence in favour of the merger hypothesis. 
However, the question if RJ is also contributing to the observed sub-mm emission remains open. These results illustrate the promising 
prospects for radio interferometry and strong gravitational lensing to study the internal structure of SMGs. 

Key words. Gravitational Lensing - Galaxies: starburst - Galaxies: high redshift - Radio continuum:galaxies - Interacting group - 
submillimetre - galaxies: individual SMM J04542 - 0301 



1. Introduction 

The detection of the cosmic infrared background (CIB) by the 
COBE satellite dPuget et all Il996t lHauser et alJ Il998h estab- 
lished that about half of the total radiation in the universe comes 
from dust -obscured galaxies th at are missing from optical sur- 
veys (see lLagache et al.l I2005L for a review). This population 
of dusty objects was first resolved by the IRAS and ISO satel- 
lites up to z ~ 1, and turned out to be dominated by lumi- 
nous and ultra-luminous IR galaxies (LIRGs and ULIRGs). The 
step into the high-z universe came with the advent of sub-mm 
and mm surveys, since the far-infrared (FIR) luminosity peak of 
high-z obscured galaxies is red-shifted into the sub-mm band 
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dFranceschini et all fl99ll: iBlain & Longahlll993l) . As a result, 
sub-mm galaxies (hereafter SMGs) turned out to be hundreds of 
times more numerous than galaxies with similar luminosities in 
the local universe, suggesting that they constitute the dominant 
contributor of the CIB and cosmic star formation at z > 1 . This 
illustrates the very important role that SMGs play in the con text 
of galaxy formation and evolution (see lLonsdale et al]|2006l for 

a review on (U)LIRGs and SMGs). 

O ne decade after their discovery with SCUBAP] (|Smail et ail 

ll997UHughes et al.lll998tlBarger et al.lfl998i: lEales et alll 19991) . 
it is generaly acepted that SMGs are heavily dust-obscured 



1 Submillimeter Common-User Bolometer Array, decommissioned 
in September 2005 from the James Clerk Maxwell Telescope on Mauna 
Kea 
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galaxies at high redshift (2 < z < 3) with ULIRG-like lu- 
minosities (Lfir ~ 10 12 L ) and star formation rates of the 
order of 1000 M Q yr~ x . This enormous bolometric luminos- 
ity seems to be dominated by star formation processes in- 
duced by galaxy interactions/mergers, althoug h a good fraction 
(~ 30 - 50%) of SMGs also host AGN activity (lAlexander et al.l 
120051) . First estimates of their physical properties indicate that 
SMGs are massive, gas rich systems (M gas ~ 10 10 - 10 11 M e ) 
in which t he starburst region has a typical scale i n the range 
1 - 8 kpc dNeri et al 1 120031; IChapman et al.l 120041: iGreve et al ' 



2005 



2008) 



iTacconi et al.l 120061: IWang et al l 120071: iBiggs & Ivisonl M 



be the progenitors of massive local ellipticals dLillv et al.l 


1999 


Smail et al.l 2002[ 20041 Webb et al. 20031: 


Genzel et al. 


2003 


Alexander et al. 2003, 2005; Swinbank et al. 


2006). 



However, these general properties of SMGs are based on 
the study of the very brightest examples of this class of object 
(S850/im <; 2 mJy), and may not be r epresentative of the e ntire 
SMG population. In fact, according to lKnudsen et al.l d2008l) . the 
dominant contribution to the sub-mm extragalactic background 
comes from the fainter (sub-)mJy sources that cannot usually 
be detected due to the confusion noise of current instruments 
(Ss50/jm ~ 2 mJy). The only way in which it has been possible 
to push bellow this sensitivity limit, is by using the lensing mag- 
nification provided by massive clusters of galaxies to increase 
the effective resolution of SCUBA. This approach has improved 
the sensitivity of sub-mm maps by factors of a few with respect 
to blank field surveys, alth ough just a handfull of faint SMGs 
have being identified so far (ISmail et al.l2002l:ICowie etal .120021: 
iKnudsenetaTl 120081) . Since so little is known about these in- 
trinsically faint sources, it is crucial to further investigate their 
properties (e.g. spectral energy distributions, morphologies, red- 
shifts, etc.) and assess whether they are different from the ob- 
served properties of brighter SMGs. 

A promissing strategy to gather information about faint 
SMGs is to study members that are multiply imaged by clusters 
of galaxies. In this cases, the magnification factor can go up to 30 
(or more), providing not only the opportunity to detect but also 
to spatially resolve the morphologies and internal dynamics of 
faint SMGs at a level of detail far great er than would otherwise 
be possible (see ISwinbank et al.l 120071 for an example of this 
technique in the optical). To date, only one multiply-imaged faint 
SMGs has been confirmed : SMM J16359+6612, located near the 
core of the cluster A2218 (iKneib et al.ll2004 120051: iGarrett et ail 
l200HlKnudsen et al.ll2008[l2009ir There are, however, two other 
clusters which seem to host multiply-imaged SMGs: A 1689 
dKnudsenet "all 120081) and MS0451.6 - 0305 dChapman et al.l 



120021: iBorvs et al.ll2004t iBerciano Alba et alj|2007l) . 

The case of MS045 1 .6-0305 is particularly interesting, since 
the "sub-mm source" (SMM J04542 - 0301) is an elongated 
(~ 1') region of 850 fim emission which is coincident with an 
optical arc and five NIR sources (see Fig. [T). While the opti- 
cal arc is the result of a strongly-lensed Lyman Break Galaxy 
(LBG) at z S pect = 2.91 1, a lens model of the cluster predicts that 
the set of NIR sources could be produced by two triply-imaged 
EROs^ located at almost the same redshift (z mo dei = 2.85 + 0.1). 
The model also indicates that the ERO pair and the LBG may 
constitute a merger at z - 2.9 in the source plane, with their 
interactio n likely being at t he origin of the observed sub-mm 
emission (IBorvs et al.ll2004l B04 hereafter). Unfortunately, the 
low resolution (~ 15" at 850 /mi) and poor positional accuracy 
(~ 2" - 3" rms) of SCUBA, makes it very difficult to confirm the 
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Fig. 1. Summary of the results presented in IBorvs et aD 
( 120041) . SCUBA 850 ^m contour map of SMM J04542 - 0301 
superimposed upon a HST image of the center of the cluster 
MS045 1 .6 - 0305. Highlighted with rectangles are an optical arc 
(ARC1) and its counter image (ARC1 ci), produced by a LBG 
at Zspect = 2.911. The red circles indicate the positions of five 
NIR sources that have been interpreted as multiple images pro- 
duced by two EROs (ERO B, lensed as B1/B2/B3, and ERO C, 
lensed as C1/C2/C3) at z mode i = 2.85 ±0.1. Assuming z = 2.9 
for both the LBG and the EROs, their predicted positions in the 
source plane would be located withing a region of ~ 10 kpc, 
suggesting that they constitute a merger. 



link between SMM J04542 - 0301 and the proposed optical/NIR 
sources. Moreover, the emission coming from the north-eastern 
and central regions of the sub-mm emission cannot be repro- 
duced by the arc and the NIR sources (see Fig. 7 in B04), sug- 
gesting that it might arise via other sources. 

A possible way to overcome this resolution problem is 
by taking advantage of the observed correlation between the 
radi o synchrotron and FIR emission in s t ar-forming galax- 
ies (Ivan der Kruillll973b ICondon et al.|[l98l iHelou et al.l ll985b 
lGarrettll2002l; lAppleton et all 12004 iBeelen et al.ll2006l). which 
also seems to hold for SMGs out to z ~ 3 ( Kovacs et al. 2006; 
IVlahakis et al1l2007l: llbar et al.ll200l; iMichalowski eTaDl2009l) . 
Thanks to this FIR-radio correlation, radio interferometric obser- 
vations can be used as a high-resolution proxy for the rest-frame 
FIR emission observed in the sub-mm. This approach has been 
extensively used to pinpoint the position of SMGs in order to 

Ivison et all 120001; 



identify their faint optical counterparts (e.g 
iBarger et alJl2000HChapman et al.ll2005b llvison et alll2007l) . 

In IBerciano Alba et alj d2007l) . we reported the detection of 
1.4 GHz radio emission coincident with SMM J04542-0301 us- 
ing VLA0 archival data. Part of this radio emission is located in 
the region between the optical arc and the ERO images, which is 
consistent with the interacting region of the hypothetical merger 
being the source of the observed radio and sub-mm emission. 
We also detected bright radio emission in the central region of 



Extremely Red Objects, photometrically defined as R - K > 5.3 



3 Very Large Array interferometer of the National Radio Astronomy 
Observatory (NRAO) 
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SMM J04542 - 0301, although it is not clear if this emission is 
produced by a high-z lensed object or AGN activity associated 
with a cluster member. 

In this paper, we present a higher resolution 1 .4 GHz radio 
map of the center of the cluster MS045 1 .6 - 0305 (MS045 1 here- 
after), obtained after combining the previous VLA archival data 
(BnA configuration) with new VLA high resolution observa- 
tions (A configuration). The details of both sets of observations 
are presented in Sect. 2, together with a description of the data 
analysis pr ocedure, which has been co nsiderably improved with 
respect to iBerciano Alba et all d2007l) . Section 3 is dedicated 
to the analysis of the final combined data set, from which we 
characterise the compact and extended radio emission detected 
within SMM J04542 - 0301 . Identification of optical (HST/ACS 
F814W) and NIR (SUBARU/CISCO K' band) counterparts for 
the radio detections, as well as their connection with the sub-mm 
emission, is discussed in Sect. 4. In Sect. 5 we describe a new 
lens model of the cluster MS0451.6 - 0305, built to investigate 
the lensed nature of the compact radio d etections. A discus sion 
about the merger scenario proposed by iBorvs et all d2004l) in- 
cluding the results from the radio observations are presented in 
Sect. 6. Summary and conclusions are presented in Sect. 7. 

The adopted cosmology corresp onds to a ACDM mo del with 
Q ra = 0.28, Q A = 0.72 and h = 73 dSpereel et al.ll2007l) . 

2. Radio Observations and Data Reduction 

The 1.4 GHz observations of the cluster MS0451 presented in 
this paper were obtained with the VLA in A and B configurations 
(A-array and B-array hereafter) using the wide-field "pseudo- 
continuum mode". This mode consists on two IFs of 25 MHz, 
centered at 1.3649 GHz and 1.4351 GHz. Each IF contains right 
and left circular polarizations and it is split into seven 3. 125 MHz 
channels. The absolute flux density scale was set by the flux cal- 
ibrator 0137 + 331 (3C48), while 0503 + 020 was used as the 
phase calibrator (PHCS hereafter). The data reduction was con- 
ducted in classic AIPS using our own semi-automatic pipeline 
written in Parselton gue (Python interface to classic AIPS, see 
iKettenis et ai1l2006l) . 

2.1. A-array Observations 

The A-array observations consist on 2 x 6 hours of data acquired 
the 5th and 10th of February 2006 (correlator integration time of 
3.3 seconds). Each IF of each observing day was independently 
self-calibrated, resulting in four data sets that were in the end 
combined to obtain the final A-array map of the target. In order 
to avoid source smearing, no averaging in frequency nor time 
was applied to the data. The different steps followed during the 
data reduction of each epoch are described next. 

First, we improved the positional accuracy of the antennas 
in our data by applying the most recent baseline corrections de- 
termined by NRAO. Corrupted data were visually inspected and 
flagged independently for each channel, IF and polarization at 
different stages of the data reduction process. 

To set the flux scale, amplitude and phase corrections were 
calculated for 3C48 using a model of this source provided by 
AIPS. Then, 3C48 was used to calibrate the bandpass shape. 
In the case of the PHCS, we first derived a preliminary cal- 
ibration in amplitude and phase using a point source model. 
Subsequent self-calibration produced a more accurate clean 
component model (CC model hereafter), which was ultimately 
used to derive the final phase calibration. After proper interpo- 
lation, the resultant amplitude and phase solutions were applied 



to the cluster data, and further refined through self-calibration in 
order to produce the optimal map for our target. 

Figure [2] (top) shows a low resolution lxl deg map of the 
cluster produced with the inner 10 kA baselines of the A-array 
data. Due to the significant number of bright sources located out- 
side the main lobe of the primary beam (black circle), 37 facets 
were used for the imaging part of the self-calibration process: 
31 overlapping facets that cover the central 0.32 deg radius re- 
gion, and 6 individual facets centered at the positions of bright 
NVSS sourcefl located in an annulus between 0.32 and 0.5 deg. 
radius. The size of each facet is 2048 x 2048 pixels, with a pixel 
scale of 0.247". Areas for cleaning were restricted by placing 
boxes around the brightest sources in each facet. The weighting 
scheme used for imaging was selected by the imagr parameter 
robust (R hereafter), with value ranging from +5 (pure natural 
weighting) to -5 (pure uniform weighting). To produce a good 
model of the compact emission, R = was used until the last 
iteration of the self-calibration procedure, where we switched to 
R = 5 to include the extended emission. 

From the sources indicated in the top panel of Fig. [2] SRC1 
is compact, and it is the brightest and closest to the target, so it 
should provide the best self-calibration solutions in the target re- 
gion. However, the extended sources SRCA and SRCB dominate 
the solutions in the short baselines. In order to get the final cali- 
bration only from SRC1 for all baselines, all the sources located 
outside the central facet were subtracted from the data using the 
final CC model of the outer 37 facets. Finally, the resultant data 
set was self-calibrated again using R = 5 during the imaging 
process. 

This procedure was followed for each IF of each epoch, re- 
sulting in four independently calibrated data sets (in which only 
the sources from the central facet remain) that were concatenated 
with dbcon to produce the final A-array map of the target. 

2.2. B-array Observations 

The B-array observations, made the 9th and 10th of June 2002, 
were retrieved from the NRAO data archive system^ (7.8 hours 
in total, correlator integration time of 10 seconds). The first 
time t hat these observations were analyzed dBerciano Alba et all 
120071) . they were treated as a single epoch and averaged on time, 
with no independent calibration for each IF. However, since our 
goal is to combine the data from both array configurations to pro- 
duce a more sensitive map of the target, the B-array observations 
were re-reduced following essentially the same procedure de- 
scribed in the previous section. The only difference corresponds 
to the self-calibration of the cluster data, which we now describe 
in detail. 

Figure |2] (bottom) shows a 1.5 x 1.5 deg low resolution map 
of the cluster produced with the inner 10 KA baselines of the 
B-array data. Note that SRC2 is brighter than SRC1, contrary 
to the situation in the A-array data (top) where SRC2 is re- 
solved. However, since SRC2 is considerably further away from 
the phase center, the best self-calibration solutions on the tar- 
get region will still be provided by SRC1. Therefore, the data 
was self-calibrated in three steps. First, a set of 7 overlapping 
facets was used to image the central 0.4 deg radius region of the 
data, while another 36 facets were centered at the positions of 
NVSS sources located in an annulus between 0.4 and 1.6 deg 



4 The NRAO VLA Sky Survey (NVSS) is a 1.4 GHz continuum 
survey covering the entire sky north of -40 deg. declination, with a 
completeness limit about 2.5 mJy 

5 Project ID AN0109, PI: Nakanishi 
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Table 1. Information about the identifyed radio sources in the core of MS0451. The columns show: position (RA, DEC), peak 
flux density (S pk), total flux density inside the 2 <x contour (S 2a), total flux density obtained by adding the flux from all the clean 
components used to model the source ( Sec), total flux density inside the 2 cr contour in the residual map (S res idX comparison 
between the total flux derived from the CC model and the 2 cr contours (^~ 2 ")> comparison between the total flux derived from 
the CC model and the gaussian fit ( 5c g~ l5cmas ) and signal-to-noise ratio (SNR). The quoted errors in RA and DEC correspond to the 
expected positional uncertainty of a point-like source in the map (FWHM/2 x SNR, were FWHM is the full with half maximum of 
the clean beam), while the error in the flux peak is estimated to be the same as the rms noise of the map. The error in the 2cr flux 
was calculated as rms x yfN, where rms is the rms noise of the map, and N is the number of beams within the area delimited by the 
2cr contour of the source. The raw "RJ bright" lists the fluxes for RJ estimated withouth taking into account the Act west extension. 
Note that the positions listed here were derived before aligning the radio map respect to the HST image. 



Name RA (+4* 54 m ) DEC (-3°) 5 Pk 5 2t r See S rcsld ^P^- s cc -s^ SNR 
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radius (big boxes in the bottom panel of Fig. [2). Each facet is 
2048 x 2048 pixels, with a pixels size of 0.814". Once the self- 
calibration procedure was completed, the CC model from the 
42 outer facets was subtracted from the data, leaving only the 
sources located in the central facet. In the second step, the central 
region of the resultant data set is imaged again using 3 1 overlap- 
ping facets of 500 x 500 pixels. After following the same self- 
calibration procedure used in the first step, the best CC model 
of the outer 30 facets was subtracted from the data. In the third 
step, this final data set is calibrated (phase only) using a solution 
interval of 5 min and R = 5. 

As in the case of the A-array observations, the four indepen- 
dently calibrated data sets obtained in this way (one per epoch 
and IF) were then combined and imaged to produce the final B- 
array map of the target. 

3. Analysis of the A+B-array map 

Before the A-array and B-array observations could be combined 
to produce the final radio map of the cluster center, it was nec- 
essary to correct for a single position shift between the two cal- 
ibrated data-sets. This correction was achieved by performing a 
global phase-selfcalibration of the B-array data set using the best 
A-array CC model. 

To produce the map with the best compromise between sen- 
sitivity and noise properties, we explored the full range of values 
available for the imagr parameter robust. The histograms of the 
resultant clean maps show that the noise distribution is closest to 
Gaussian when R = 0.6 (see Fig - [3j ■ Note also that the rms noise 
of the R = 0.6 map (~ 10 fjjy) is only 2% higher than the one 
obtained for R = 5, which is the most sensitive map that can be 
produced with this data. 



3.1. Catalog of radio sources 

Figure|4]shows a detail on the R = 0.6 map of the cluster center 
produced by cleaning the data down to a depth of 5.5 fi}y. The 
contours indicate the region of extended sub-mm emission re- 
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Fig. 3. Histogram of the A+B array clean map produced using 
R = 0.6. The black curve is the best Gaussian fit to the noise 
peak, while the positive excess on the right side is due to source 
emission. The good agreement between the fit and the data in the 
negative side of the gaussian indicates that the noise of the map 
is well behaved, making faint detections more reliable. 



ported by B04, in which six of the eight identified radio sources 
are located. Postage stamps of these six sources and their CC 
models are shown in Fig. [5] 

The radio detections located outside the sub-mm emis- 
sion, labeled Fd and Fe, are the same sources reported in 
iBerciano Alba et al.l d2007l) . Their morphology in the new high 
resolution radio map is very elongated, typical of radio jets pro- 
duced by AGNs. In the case of the other 6 sources, the mor- 
phology and distribution of their clean components show that, 
with the exception of the (almost) point-like structure of El, the 
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Table 2. Gaussian fitting results for the identified compact radio sources in the core of MS0451. The columns show: position 
(RA, DEC), peak flux density (Spk), total flux density (S Gauss), deconvolved Gaussian sizes (major axis, minor axis and position 
angle) with their corresponding formal errors, and the sum of the square of the residuals reported from the fit (SSR). Cases where 
the parameters of the Gaussian fits are not well constrained are indicated by a dash. Note that the positions listed here were derived 
before aligning the radio map respect to the HST image. 
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70.5 ± 39.0 


3.4 
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Fig. 6. Detail of Fig.|4] The grey scale has been modified to make 
more evident the low level background in this region of the radio 
map. 

rest of the sources are extended. In particular, RJ shows a bright 
elongated main body and an irregular 4<x extension towards the 
west. Note however, that the amorphous extensions present in 
the 2<x contours of El, CR2 and E3 are not considered to be part 
of their real structure. In the following, references to the 2cr con- 
tours of these sources will implicitly exclude these amorphous 
extensions. 

The peak positions and peak fluxes of all the radio detections 
were calculated with the task maxfit, which fits a quadratic func- 
tion to a rectangular area that encloses the source (in our case, a 
rectangle delimited by the 2<x contours of the source). The total 
flux, on the other hand, was determined in two different ways: 
(i) by integrating the emission over the area defined by the 2<x 
contour of the source (S 2a-), and (ii) by adding the flux of all 
the clean components of the source model (Sec)- In this way, we 
can calculate the flux excess contained in the clean component 
model respect to the 2<x flux, which is an estimate of the amount 
of flux contained in the wings of the PSF-convolved sources be- 
low the 2<t contour. 

The results obtained for the different parameters are listed 
in Table Q] Note that, since the cleaning procedure is stopped at 
certain flux level (half a sigma below the noise level), it is ex- 



pected to have some left-over emission from the sources in the 
residual map, specially in the case of extended sources. To esti- 
mate this left-over flux, we have integrated the emission inside 
the 2<t contours of each source in the residual map (see S les id in 
TableQ}. The results show that the S cc for CR1 and RJ is under- 
estimated by ~ 6%, and by ~ 10% in the case of E2 and E3. The 
negative values for 5 res id obtained for El and CR2 indicate that 
all the flux has been included in the CC model. 

We also tried to characterize the properties of the six sources 
located within the sub-mm emission region by comparing them 
with Gaussians. The task jmfit was used to fit simultaneously all 
the parameters of an elliptical Gaussian in the same rectangular 
area used with maxfit. Since the shape of RJ is far from being 
elliptical , we restricted the Gaussian fitting to the area that only 
contains the 2<x contours of the bright main body. The center and 
FWHM of the best fitting gaussians are shown in Fig. [5] 

The parameters derived from the best fits and their formal 
errors are listed in Table [2] where the sum of the square of the 
residuals is a measure of the relative goodness of the fit. As ex- 
pected, the worst fits are obtained for RJ and E3, since they have 
the most irregular shapes. The best fits correspond to CR1 and 
E2, but note that the FWHM of the Gaussian in CR1 does not 
include the north-east 2<x extension. 

When compared with Table [TJ we see that the Gaussian 
fit provides basically the same estimates for the position (RA, 
DEC) and flux peak (Spk)- The nominal integrated flux (S Gauss) 
is in general larger than Sec (except for El), but consistent 
within the errors. The largest discrepancies are found for E3 and 
CR2 (see Scc s ^ SGm " i n Table [T), probably because S Gauss is in- 
cluding some flux from the 2<x amorphous extensions in these 
two sources that were not accounted for in the clean component 
model. 

Finally, we would like to emphasize that, despite the faint- 
ness of CR1 and CR2 (SNR ~ 4.6), the histogram presented in 
Fig. [3] and the distribution of -4<x peaks in the final map (only a 
small number, located far away from the target region) suggests 
that they are real detections. 

3.2. Disentangling the compact and extended radio emission 

As can be seen in Fig. [6] the region between El and E2 shows 
an enhanced low level background when compared with the rest 
of the map, suggesting the presence of diffuse radio emission 
that has been partially resolved out. From now on, we will re- 
fer to this diffuse emission as the extended component, whereas 
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El, E2, E3, CR1 and CR2 will be consider compact sources by 
comparison. 

If present, this extended component would only be detectable 
on short baselines, while the emission produced by compact 
sources will be picked up by all baselines at different resolution 
levels. Therefore, a map produced using only the short baselines 
does not allow us to distinguish between extended emission and 
a compact source observed at low resolution. The only way to 
isolate the extended component is by subtracting a model of the 
compact sources from the data. 

Given the faintness of CR1 and CR2, and the complicated 
low level background in this region, the choice of boxes during 
the cleaning procedure becomes a highly subjective issue, and 
even the decision of cleaning CR1 and CR2 (and hence their CC 
models) can be questionable. For this reason, we produced a new 
version of the map presented in Fig. [4] (following the same pro- 
cedure) in which CR1 and CR2 have not been cleaned. Then, we 
calculated Sp^, S n a and Sec for all the sources in the new map, 
and compared them with the values presented in Table [T] In the 
case of E2 (which is the closest source to CR1 and CR2) we find 
differences of ~ 10% for Sec and ~ 14% for S20-, whereas for 
the rest of the sources the differences in all three parameters are 
less than 7%. Since these variations are well withing the errors 
estimated for Spk and S n a in Table[T] we can use the CC models 
of CR1 and CR2 presented in FigfS] without compromising the 
results. 

Fig. [7] shows three ~ 5" resolution maps of the region be- 
tween El and E2, produced by tapering the data with a Gaussian 
of FHWM = 30 k/l. In map (A) no compact sources were sub- 
tracted before the tapering and imaging. The central map (B) 
shows the result of subtracting the CC model of El and E2 from 
the data, and map (C) shows what is left after subtracting all 
the compact sources. Given the robust detection of a 6.6cr elon- 
gated source in panel (B), the presence of resolved extended 
radio emission in this region of the ~ 2" map (Fig. [6j is con- 
firmed. Using the 2<x contour of map (B) as template, the inte- 
grated fluxes in panels (B) and (C) are 5b = 173.52 + pjy and 
Sc = 84.14 ± piy. The flux error was calculated as rms x Vw, 
where rms is the rms noise of the map, and N is the number 
of beams within the area delimited by the 2a contour of map 
(B). These numbers indicate that ~ 50% of the flux in panel 
(B) has been included in the ~ 2" CC model of CR1 and CR2. 
Therefore, it seems that CR1 and CR2 constitute two relatively 
compact regions of an extended radio source located between 
El and E2. However, given the faintness of CR1, CR2 and the 
extended component in panel (C), deeper observations will be 
required to confirm this result. 

Finally, to estimate the relative contribution of all the com- 
pact sources respect to the extended emission, we used the 2cr 
contour of map (A) as template to calculate the integrated flux 
of each map in Fig. [7] S A = 285 + 33 pjy, S B = 204 ± 33 pjy 
and Sc = 109 + 32 piy. The comparison of these fluxes indicates 
that 38% of the emission comes from an extended component, 
whereas 62% is produced by the compact sources (28% from 
El and E2, 34% from CR1 and CR2). As it will be discussed in 
detail in Sect. 14.11 the presence of this extended radio compo- 
nent constitutes a strong observational evidence in favor of the 
merger scenario proposed in B04. 

4. Multiwavelength Counterparts 

Table [3] summarizes the relevant information about all the mul- 
tiwavelength data of MS0451 (optical, NIR and sub-mm) that 
we have collected from the literature to compare with our radio 



map. To identify possible counterparts of the radio sources, the 
radio and K' band observations were aligned with respect to the 
HST imag^B using the iraf package ccmap. For the alignment of 
the radio map we used 13 compact radio sources with reliable 
bright counterparts in the HST image, whereas 93 compact NIR 
sources were used to align the K' band image. We find an rms 
scatter of 0.29" between the radio map and the HST image, and 
0.06" between the K' band and the HST images. Therefore, the 
inferred error in the alignment is 0.29/ Vl3 = 0.08" for the ra- 
dio map and 0.06/ V93 = 0.006" for the K' band image. These 
errors are negligible compared with the (relative) astrometric ac- 
curacy of the sources in each map, so they will be ignored in 
further analysis. 

4. 1 . Comparison between radio and optical/NIR 

Figures [8] [9] and [14] show different details of the HST image 
where the radio detections are located. The positions of the rele- 
vant sources are indicated with circles: blue for the radio sources 
and red for the K' band sources. The counterparts of the sub- 
mm emission proposed by B04 have been indicated following 
the nomenclature introduced in Fig.Q] Note that the images B3 
and C3 are supposed to be unresolved in the K' band image pre- 
sented in B04, and therefore referred to as B3/C3. 

The relative position of the radio and NIR sources in these 
figures suggests that El, E2 and E3 might be associated with 
Bl, B2 and B3/C3. To better quantify these associations, as well 
as identify other possible counterparts, we have used Bayesian 
inference to calculate the probability that the compact radio 
sources are physically associated with their nearest optical/NIR 
sources. The mathematical expression of this probability has 
been derived as follows. 

Lets P(p\d) be the probability density of the true distance p 
between two sources given an observed distance d. Following 
Bayesian inference, this probability is given by: 



P(fl\d) 



urn m 

P(d) 



(1) 



where L{d\p) is the likelihood of the observed distance given 
a true distance, P(p) is the probability density of the true distance 
before any measurement, and P(d) is the probability density of 
the observed distance for all possible distances (a normalization 
constant). 

As pointed out in lChurchman et al. (2006), the likelihood of 
the observed distance between two sources which positions are 
described by a Gaussian distribution is not Gaussian. Instead, 
this likelihood in two dimensions is given by: 



r/j| % dkidfi/io-l + crl)) 
L(d\p) = exp 



oi + oi 



-p 2 -d 2 
{2(o- 2 + o- 2 2 ) 



(2) 



Where cr\ and <x 2 are the errors in the position of each source, 
and Io is the modified Bessel function of integer order zero. 

Since the area inside an annulus at radius p increases as p 
dp, we assume that, for a source located at a random position, 
P(p) = p. Therefore: 



P(p\d) cc L{d\p) ■ p 
Which normalized gives: 



B{p\d) = 



L(d\p)-p 

r-oo 

Jo L(d\p)dp 



(3) 



(4) 



6 astrometry matched to S ubaru/Suprime-Cam observations of 
MS0451 for the PISCES survey (iKodama et all2005l) 
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Table 3. Observations details. 



Instrument 


Frequency (filter) 
Gum) 


Observing dates 
(dd/mm/yy) 


Total exposure time 
(ks) 


depth 
(magnitudes) 


resolution 

(arcsecs) 


references 


HST/ACS 


0.83 (F814W) 




2.4-2.7 per pointing 


25 


1.0 


Moran et al. (2007) 1 


SUBARU/CISCO 


2.13 (K') 


17-18/11/2000 


9.22 


22.1 


0.6-0.8 


Smith et al., in prep. 
Takata et al. (2003) 


JCMT/SCUBA 


850 


03/09/1998 


22 




15.0 


Borvs et al. (2004) 



1 Data products from their wide field survey in MS0451 are available at http://www.astro.caltech.edu/ smm/clusters/ 
Further information on the data and survey can be found in Chapter 2 of the PhD thesis of Sean Moran. 



Table 4. Conunterparts of the radio detections. The columns 
show: possible counterparts of the radio detections (Optical/IR 
source), mesured offset between the radio source and the op- 
tical/IR source (af bsX radius of the optical/IR source (r src ) and 
probability that the centroid of the radio source lies inside the 
area of the optical/IR source (P assoc .). The coordinates of the 
counterpart sources can be found in tables [5] and [6] Sources 
which positions were taken from the K' band image are indi- 
cated with an asterisk (otherwise we used the HST coordinates). 
The positions of the radio sources used to calculate c/ b s corre- 
spond to the values derived with maxfit in the aligned radio map. 



Radio 


Optical/IR 


dobs 


r src 


p 

1 assoc. 


source 


source 


(arcsecs) 


(arcsecs) 


(%) 


El 


Bl* 


0.38 


0.6 


67.8 


CR1 


ARC 1.1 


1.12 


0.15 


< 1 




ARC 2.1 


1.40 


0.25 


< 1 




ARC 3.1 


1.60 


0.20 


< 1 


CR2 


G2 


1.39 


0.70 


< 1 




C2* 


1.82 


0.60 


< 1 




ARC 1.2 


1.92 


0.15 


< 1 


E2 


B2* 


0.26 


0.60 


79.29 




G10 


0.74 


0.23 


< 1 


E3 


4.3* 


0.77 


0.60 


20.7 




5.3* 


1.22 


0.30 


< 1 




G23 


1.54 


0.35 


< 1 


RJ 


G8 


1.91 


0.30 


< 1 




Gil 


2.40 


0.26 


< 1 




G14 


1.39 


0.60 


< 1 



The integral of this function between fii and /i2 gives the 
probability that, for a given pair of source centroids with ob- 
served offset d, the true distance between them lies between 
and jU2- Intuitively, one would say that the condition for two 
sources to be counterparts of each other is fi — 0. The problem 
is that, in the case of extended sources like galaxies, the peak 
of emission at different wavelenghts can be located at different 
positions, which means that [i + despite the fact that they are 
physical counterparts. For this reason, we will consider that a 
radio source has an optical/NIR counterpart when the centroid 
of the source emission lies inside the area of the optical/NIR 
source. The probability of meeting this condition is given by: 

1 assoc. — BiiAd) dn (5) 

Jo 

A list with the nearest optical/NIR sources of each radio de- 
tection, together with their associated probabilities, is presented 
in Table [4] As expected, only the images of ERO B have a non- 
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Fig. 9. Detail of the HST image in the region of RJ, including 
the critical curve at z = 2.911 predicted by the lens model de- 
scribe in sect |4.1| (black line). The contours correspond to the 
~ 2" resolution radio map presented in fig|4] The radius of the 
circles show the estimated positional uncertainties: 0.1" for RJ 
(see Table [Qi and ~ 0.08" (half a pixel) for G14. The blue line 
indicates the approximate orientation of the major axis of the 
radio source. For comparison, the yellow line shows the orien- 
tation of the RJ centroid respect to the galaxy G14. The black 
ellipse indicates the beam size of the radio map (2.14" x 1.71" 
at a position angle of - 1 .43). 

negligible probability of being associated with some of the radio 
sources (El, E2 and E3). However, if El, E2 and E3 are indeed 
the radio counterparts of Bl, B2 and B3, they have to be multi- 
ple images produced by a radio source at z ~ 2.9. A quantitative 
analysis of this lensing scenario is presented in Sect. 15.41 

Finally, note that no optical/NIR counterpart has been identi- 
fied with RJ. Therefore, given the depths of the optical and NIR 
images, it seems more likely that the extended morphology of 
RJ corresponds to an AGN rather than a resolved low redshift 
star forming galaxy. Within the AGN scenario, we expect that 
the peak of the radio emission corresponds to the position of an 
undetected optical source, and the extensions in opposite direc- 
tions are two jets coming from it. Another conceivable scenario 
would be that the galaxy G14 is an AGN host with a one sided 
radio extension. However, the mayor axes of RJ does not seem to 
be aligned with G14, and this kind of AGN is not very common. 

4.2. Comparison between radio and sub-mm 

In this section, we will use the FIR-radio correlation to check 
whether the radio detections could be associated with the ob- 
served sub-mm emission. To make this kind of analysis, it is 



8 



Berciano Alba et al.: Internal radio structure of a faint sub-mm galaxy 




23456789 



Fig. 11. Change of the parameter with redshift for different 
characteristic dust temperatures. The horizontal dashed lines in- 
dicate the minimu m and maximum valu e of found in the SMG 
sample studied bv lKovacs et al.l (|2006). The dotted vertical line 
indicates the position of z-2.9. 



necessary to match the resolutions of the radio and sub-mm ob- 
servations. Therefore, the A+B array data was tapered with a 
Gaussian of 8.5 and restored with a clean beam of 15" x 15" 
at the end of the imaging process. The left panel of Fig.fTOlshows 
the final tapered radio map (white contours) superimposed upon 
the sub-mm map (black contours) and the HST image of the 
cluster core. The positions of the detections in the high reso- 
lution radio map have been indicated with crosses as reference. 

If all the observed sub-mm emission would be produced by 
a single SMG that does not host an AGN, it is expected that 
the radio and sub-mm morphologies would resemble each other. 
This is because the origin of the FIR-radio correlation seems to 
be linked to massive star formatiotfl which means that both the 
radio and sub-mm emission are originated in (approximately) the 
same regions of the galaxy. However, the 5 850 ^ m /S 1.4 GHz 

flux 

ratio observed in sub-mm galaxies displays a broad scatter which 
strongly de pends on the characteristic dust tempera ture and the 
redshift (e.g lBlain et a"T]|2002l: IChapman et alj|2005l) . Therefore, 
if the observed sub-mm emission is being produced by several 
blended SMGs, we might find "morphological inconsistencies" 
like the one present in the left panel of Fig. [10] (note that the 
brightest peak of the radio emission is located in the region of 
RJ instead of been coincident with the brightest sub-mm peak). 

Interestingly, if the CC model of RJ is subtracted from the 
data, the morphology of the radio emission becomes remarkably 
similar to the sub-mm emission (right panel of Fig. [TOt . This 
result strongly suggests that the source/sources responsible for 
the radio emission observed in this panel are also responsible for 
the bulk of the sub-mm emission. The other conclusion derived 
from this comparison is that, if RJ is contributing to the sub-mm 
emission, either its properties (redshift and/or dust temperature) 
are different from the properties of the source/sources associated 
with the other radio detections, or it has a "radio excess" due to 
an AGN. 

Using 350 /mi, 850 //m and 1.4 GH z observations of 15 
bright SMGs with spectroscopic redshifts, K ovacs et alj d2006l) 

7 which corresponds to a clean beam of 15.78" x 14.65" 

8 while young massive stars produce UV radiation that is re-emitted 
in the FIR by the surrounding dust, old massive stars die as SN produc- 
ing electrons that are accelerated by the galactic magnetic field gener- 
ating synchrotron radio emission. 



(K06 hereafter) conclude that the FIR-radio correlation remains 
valid for SMGs with z ~ 1 - 3 and luminosities between 
10" - 10 13 M Q (except when they are radio-loud AGN). To 
see how our system compares with this sample, we have made a 
more quantitative analysis of the FIR-radio correlation using the 
q^ parameter introduced in K06: 

qL = i og ( ) (6) 

qh *\[4.52THz]Li.4ghz/ 

were total FIR luminosity is inferred from the flux density at 
850 /on assuming a grey body model: 



L FIR = 4nD 2 L r(4- +0)a4 +P) {e hY/m - l) vS 85 o A « m (7) 

and the radio luminosity can be derived as: 

LiAom = 4tt£>2 s ughz (1 + z)^ 1 (8) 

In these formulas, Z\ is the luminosity distance, /3 is the ef- 
fective emissivity index (we assume /3 = 1.5), is the charac- 
teristic dust temperature, v corresponds to the observed sub-mm 
frequency, and a is the radio spectral index (we assume a = 0.7). 

At the position of the sub-mm peak, the radio and sub-mm 
fluxes were derived using the area delimited by the 3<x radio con- 
tour shown in the left panel of Fig.[10](5'850//m = 15.5 + 2.8 mJy 
and S 1.4GHz = 241 + 42 /dy). The different values of q^ de- 
rived from these fluxes as function of redshift and I'd are shown 
in Fig. QT| The horizontal dashed lines indicate the minimum 
and maximum value of q^ found in the K06 sample (discard- 
ing AGN hosts). Note that there is a wide range of z and T& for 
which the observed fluxes are consistent with K06. If we now 
assume that the radio and sub-mm emission are produced by a 
(lensed) star forming galaxy at z = 2.9, only temperatures be- 
tween ~ 20 - 40 K would be allowed, which is the temperature 
range in which most of the K06 sources lie. 

Since AGNs do not follow the FIR-radio correlation, an esti- 
mate of the <7l value for RJ could provide extra evidence to con- 
firm/discard the possible AGN nature of this source. However, 
that would require a higher resolution 850 /mi map, to confirm 
the connection of RJ with a discrete sub-mm source and get an 
accurate estimate of its sub-mm flux. In addition, the redshift 
of RJ needs to be determined in order to break the degeneracy 
between z and Td illustrated in Fig.fTTI 

Finally, is important to mention that, although the energy 
output of the majority of SM Gs is dominated by sta r formation, 
about 30 - 50% host an AGN ( I Alexander et al.l2005l) . Therefore, 
even if the source associated with RJ is an AGN, it could still 
be contributing to the sub-mm emission observed in that region. 
Other potential contributors to the observed sub-mm emission 
are indicated in the right panel of Fig.fTOl 

5. Gravitational lens modeling of MS0451 

In order to investigate the possible lensed nature of the radio de- 
tections, we used the publicly available LENSTOoJ^code to model 
the mass distribution of the cluster. This modeling involves an 
optimization procedure, aimed to find the mass model parameter 
values that best reproduce the observational constrains (positions 
of the multiply imaged systems). For a detailed description of 

9 http://www.oamp.fr/cosmology/lenstool/ 
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Table 5. Galaxy catalog for the new lens model of MS0451. The colums show: galaxy number (ID) galaxy coordinates (RA, 
DEC), semi-major axes (a), semi-minor axes (b), position angle (9) and magnitude in the HS T image (mag). An asterisk in the ID 
number indicate galaxies that were not included in the model presented in lBorvs et al.l (120041) . 
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(arcsec) 


(arcsec) 
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/"PS 1 AW7\ 
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01* 


12.933 


00 


41.494 


0.24 


0.23 


65 


21.9511 


02 


12.884 


01 


01.056 


0.70 


0.55 


40 


20.1111 


03* 


12.886 


00 


55.685 


0.12 


0.09 


-28 


24.6286 


04 


12.855 


01 


09.205 


1.50 


0.60 


59 


19.0257 


05 


12.831 


01 


03.476 


0.17 


0.15 


55 


23.2229 


06* 


12.723 


00 


42.529 


0.24 


0.17 


30 


22.3768 


07 


12.697 


00 


54.889 


0.53 


0.40 


140 


20.8427 


08 : 


12.655 


01 


16.586 


0.30 


0.10 


-44 


22.5806 


09 


12.639 


01 


04.288 


0.70 


0.25 


-80 


20.7576 


10* 


12.634 


01 


00.949 


0.23 


0.17 


50 


22.4744 


u* 


12.610 


01 


: 15.312 


0.26 


0.16 


-69 


22.1831 


12* 


12.563 


01 


04.426 


0.19 


0.19 


20 


22.5221 


13 


12.492 


01 


08.202 


0.90 


0.50 


81 


19.6345 


14 


12.461 


01 


17.075 


0.60 


0.45 


105 


20.0909 


15* 


12.334 


01 


15.892 


0.25 


0.25 


-15 


21.9945 


16 


12.283 


00 


53.001 


0.55 


0.45 


50 


20.3429 


17 


12.280 


00 


: 54.816 


0.64 


0.42 


115 


20.0471 


18 


12.189 


01 


: 01.348 


0.75 


0.73 


10 


19.5202 


19 


11.962 


01 


: 04.756 


0.70 


0.46 


63 


20.1070 


20 : 


11.833 


00 


: 59.879 


0.45 


0.25 


16 


21.1736 


21 


11.817 


00 


: 57.484 


0.53 


0.45 


15 


20.3933 


22 


11.312 


01 


: 01.400 


0.55 


0.45 


66 


20.6069 


23 


11.032 


01 


: 26.400 


0.35 


0.30 


-40 


21.4960 


24 


10.884 


00 


: 57.231 


1.20 


0.85 


50 


18.1360 


25 


10.853 


00 


: 51.865 


1.30 


1.00 


30 


18.1454 


26* 


10.767 


01 


: 24.772 


0.12 


0.08 


21 


25.2539 


27 


10.076 


00 


: 50.338 


0.47 


0.47 


56 


20.5351 



the strong lensin g methodology f ollow e d in the lenstool soft- 
ware we refer to Limousin et al.l (120071) . lJullo et all (120071) and 
the appendix A2 of ISmith et al.1 (120051) . 

The current best lens model of the cluster MS0451, pub- 
lished in B04 (B04 m odel hereafter), wa s produced with the first 
version of lenstool (iKneib et al.lll993h . This version is based 
on a downhill x 2 minimization (which can be very sensitive to 
local minimum in the likelihood distribution), and does not pro- 
vide estimates of the er rors on the optim ized parameters. The 
latest lenstool version (lJullo et alj|2007l) used in this paper in- 
cludes a Bayesian Monte Carlo Markov chain optimization rou- 
tine, which allows to determine errors in the optimized parame- 
ters and lowers the probability of ending in local x 1 minimum. 

5. 1 . Mass Model 

The cluster mass model is constructed using two different com- 
ponents: (i) a cluster halo, which represents the DM component 
on cluster scales and the baryonic intracluster gas, and (ii) 27 
galaxy halos that account for individual galaxie s. Both compo- 
nents are described using a P IEMD profile (e.g. iLimousin et al.l 
l2005t lElfa sdottir et al. 2007), parametrized by position (RA, 
DEC), position angle (ff), ellipticity (e), core radius (R CO i e ), scale 
radius (R cut ) and velocity dispersion (cr () ). The density distribu- 
tion for this profile is given by : 

p(r) = 8° (9) 

The galaxy-scale component contains the cluster members 
located relatively close to the area were the multiple images 



are formed, because they are the ones that have the strongest 
effect in the potential of that region. We also included some 
fainter galaxies located close to the strong lensing constraints, 
since they are known to perturb the strong lensing configuration 
dMeneghetti et al.ll2007l) . The details of our galaxy catalog are 
listed in Table [5] and their positions and sizes are displayed in 
Fig. [12] Note that the galaxy catalog used in the B04 model (39 
galaxies in total) includes cluster members located farther away, 
but not the faint galaxies mentioned before. 

5.2. Multiple images 

Figure Qj] shows a close up of ARC1 as seeing on the ACS im- 
age. From its knotted structure we can identify 3 sets of mir- 
ror images, each of which should have its own counter-image. 
However, since it is not possible to distinguish the different knots 
in the de-magnified image of the arc (ARC1 ci, see Fig.fT4li. we 
assigned the same position to the counter-image of each set of 
constraints coming from the arc. 

The positions for the EROs use d in our model com e from the 
Subaru K' band image published in lTakata et al.l d2003h . which is 
deeper and has higher resolution than the CFHT K' band image 
published by B04. Based on their model, B04 identified the NIR 
source labeled B3/C3 as the result of the blending between the 
counter-images of B and C. However, we have detected a faint 
source (labeled '5.3' in Fig. [T4l) located very close to B3/C3 
in the Subar u K' b and image. This source was not reported in 
iTakata et al.l d2003l) because it is too faint to be classified as a 
(Tadafumi Takata, private communication). Since it is 

10 Distant Red Galaxies, photometrically defined as J — Ks > 2.3 
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Fig. 13. Detail of the optical arc as seeing in the HST image. The 
circles indicate the positions of the three sets of mirror images 
listed in Table [6] 

Table 6. Constraints for the new lens model of MS0451. 

Sources used as c onstraints in the lens model presented by 
iBorvs et all d2004l) are indicated with an asterisk in their ID 
number. 



Source name 


ID 


RA (+4'' 54"') 


DEC (-3°) 


Band 






J2000 (sec) 


J2000 (',") 






1.1* 


12.956 





53.379 


F814W 


ARC1 


1.2* 


12.931 





58.029 






1.3* 


11.124 


1 


26.641 






2.1 


12.958 





54.869 


F814W 


ARC1 


2.2 


12.938 





57.446 






2.3 


11.124 


1 


26.641 






3.1 


12.958 





55.147 


F814W 


ARC1 


3.2 


12.945 





56.710 






3.3 


11.124 


1 


26.641 




Bl 


4.1* 


12.806 





44.344 


K' band 


B2 


4.2* 


12.684 


1 


01.488 




B3 


4.3* 


10.927 


1 


24.795 




CI 


5.1 


12.822 





47.834 


K' band 


C2 


5.2 


12.747 





59.481 




C3 


5.3 


10.897 


1 


25.627 




RC1 


6.1 


12.890 





53.907 


1.4 GHz 


RC2 


6.2 


12.868 





59.696 





possible that this source is the resolved counter-image of ERO 
C, we decided to include it as a constraint in the new model. The 
full set of constraints used for the optimization is listed in Table 

m 
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Fig. 14. Detail of the source E3, including its radio contours as 
presented in figE] Top panel: HST image of the region, in- 
cluding the positions of the counter images of the optical arc 
(ARC ! ci) and the ERO pair (B3/C3) reported by lBorys et all 
d2004l) . Bottom panel: K' band image of the region, indicat- 
ing the constraints used in the lens model (4.3 and 5.3). The 
radius of the circles show the estimated positional uncertainties: 
~ 0.20" for the K' band sources (astro metric error for the stan- 
dard stars quoted in lTakata et al. 1 120031) . 0.2" for E3 (see Table 
[TJ and ~ 0.08" (half a pixel) for the optical sources. The black 
ellipse indicates the beam size of the radio map (2.14" x 1.71" 
at a position angle of — 1°.43). 



Table 7. Ranges in which the free parameters are allowed to vary 
during the lens model optimization. 



Parameter 


Opt. range 


X 


±6" from BCG center 


Y 


±6" from BCG center 


e 


0-0.7 


e 


0- 180 deg 


00 


900- 1800 km/s 


^core 


30- 120 kpc 


°° 


120- 140 km/s 




30- 180 kpc 



5.3. Optimization 

The set of free parameters used during the optimization are: (i) 
all the parameters that characterize the cluster halo (except /? C oi- e ) 



and (ii) the velocity dispersion (cr*) and scale radius (R* ut ) of a 
galaxy at z = z c i uster with a typical luminosity L* that corre- 
sponds to an observed magnitude of m = 18.8. 
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Table 8. Most likely mass model parameters. The columns show: position respect to the brightest cluster galaxy (X,Y), ellipticity 
of the mass distribution (e, expressed as a 2 - b 2 /a 2 + b 2 , were a and b are the semi-major and semi minor axes of the ellipse that 
describes the light distribution of the clump), position angle (8), velocity dispersion (cr ), core radius (r C0le ) and scale radius (r cut ). 
Error bars correspond to lcr confidence level as inferred from the mcmc optimization. Values into brackets are not optimized. When 
the posterior distribution is not Gaussian, we report the mode and asymmetric error bars. The coordinates of the brightest cluster 
galaxy are RA = 4 h 54 m 10.87 s , DEC = -3°0'54.00". 



Clump 


X 

(") 


Y e 

(") 


e 

(deg) 


(To 


'< core 

(kpc) 


''cut 

(kpc) 


Cluster 


-0 5 +4J 


2.7+1.5 0.69!™ 


31.3±1.1 


1144.8±36 


60.8+24 


[1 500] 


L* elliptical galaxy 








135.1+43.9 


[0] 


3o.r« 9 6 



Table 9. Best mass model. The columns show: position respect to the brightest cluster galaxy (X,Y), ellipticity of the mass distri- 
bution (e, expressed as a 2 - b 2 /a 2 + b 2 , were a and b are the semi-major and semi minor axes of the ellipse that describes the light 
distribution of the clump), position angle (ff), velocity dispersion ((To), core radius (r core ) and scale radius (r cut ). Values into brackets 
are not optimized. The coordinates of the brightest cluster galaxy are RA = 4 h 54 m 10.87 s , DEC = -3°0'54.00". 



Clump 


X 


Y 


6 e 


O"0 


'"core 


r cut 




(") 


(") 


(deg) 


(km/s) 


(kpc) 


(kpc) 


Cluster 


-0.771414 


2.256631 


0.695267 31.146240 


1117.322994 


7.667450 


[1 500] 


L* elliptical galaxy 








150.421411 


[0] 


3.880257 



The ctq and R cut of all the galaxy halos included in the mass 
model are derived from the luminosity of their associated galaxy 
using the following empirical scaling relations: 



(L\ l/2 I L\ 



1/4 



(10) 



The scaling relation for ctq assumes that mass traces light, 
and its origin resides in the Tully-Fisher and Faber-Jackson rela- 
tions. The scaling relation for the radial parameter assumes that 
the mass-to-light ratio is constant for all galaxies. The zero point 
of these relations is set by the magnitude of the L* galaxy. The 
rest of the parameters of all galaxy halos (RA, DEC, 9 and e) are 
fixed to the values measured from the light distribution of the 
related galaxy, while their R COK is set to zero. 

The ranges in which the free parameters are allowed to vary 
during the optimization are listed in Table [7] Note that the scale 
radius describes the properties of the mass distribution on scales 
much larger than the radius over which the multiple images can 
be found. Therefore, since strong lensing cannot give any re- 
liable constrains on this parameter, its value is set to an arbi- 
trary large number (1500 kpc). The ranges adopted for the L* 
gala xy are motivated by galaxy-ga l axy lensing studies in clus- 
ters dNataraian et al.lll998ll2002ara : lLimousin et al.ll2007l) . 

Since the redshift of the EROs has not been confirmed spec- 
troscopically, the modeling procedure was carried out in four 
steps. In the first step, only the sets of images from ARC1 at 
Zspect = 2.91 1 where used as constraints in the optimization. The 
resultant best model was then re-optimized, including the con- 
straints provided by ERO B but leaving Zb as free parameter. As 
a result, the predicted redshift of ERO B is z B = 2.93 + 0.13. In 
the third step, the constraints from ERO C were used in a new 
re-optimization of the best model obtained in the previous step, 
assuming Zb = 2.9 and leaving zc as free parameter. This pro- 
vides a zc = 2.85 + 0.06, which is the same redshift that B04 de- 
rived for the ERO pair. In the final step, the previous best model 



Table 10. Results of the best models for each set of constraints. 
The columns show: ID of each set of constraints (System), re- 
duced x 2 obtained for each model (x 2 final, zfix and zfree) and 
average separation between the observed and predicted position 
of the lensed images for each model (rms final, zfit and zfree). 



System 


2 

X 


2 

X 


x 2 


rms 


rms 


rms 




MFINAL 


MZFIX 


MZFREE 


MFINAL 


MZFIX 


MZFREE 


1 


3.89 


6.71 


4.04 


0.34 


0.45 


0.35 


2 


2.02 


0.68 


1.28 


0.25 


0.14 


0.20 


3 


1.33 


0.74 


0.95 


0.20 


0.15 


0.17 


4 


3.86 


8.53 


6.32 


0.34 


0.51 


0.44 


5 


0.50 


1.48 


0.62 


0.12 


0.21 


0.14 


6 




35.46 


3.70 




1.26 


0.41 



is re-optimized (in the image plane) including all the constraints 
wifhz = 2.911. 

The Bayesian MCMC optimization routine included in the 
new version of lenstool provides two kinds of outputs: (i) 
the likelihood of reproducing the observed constraints, indepen- 
dently derived for each free parameter of the mass model, and (ii) 
the set of model parameters that provides the best fit to the input 
data. The most likely values for the model parameters obtained 
from these histograms are listed in Table [8j whereas the param- 
eters of the best model (mfinal hereafter) are listed in Table [9] 

The image positions are well reproduced, with image plane 
positional rms differences between ~ 0.1" and 0.3" (see Table 
[Tol l. The projected mass within the Einstein radius (here approx- 
imated by the ARC1 distance from the center, is M2D (30") = 
1.7310^0. 



5.4. Analysis of the radio data using the new lens model 

Since CR1 and CR2 are located in a region with several 
multiply-imaged systems (see Fig. O, it is conceivable that they 
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Table 11. General results of the best models. The columns show: 
Nickname of the model (Model), sum of the reduced^ 2 of all the 
optical/NIR sets of constraints (x 2 optical), reduced x 2 for the 
radio set of constraints (x 2 radio), average separation between 
the observed and predicted position of all the lensed images used 
as constraints (aver rms), probability of the model to reproduce 
the observations (log(evidence)). 



Model x 1 optical x 1 radio aver rms log(evidence) 

mfinal 11.60 0.26 -24.503 

mzfix 18.14 35.46 0.53 -45.340 

mzfree 13.21 3.70 0.30 -29.741 



The results (shown in Fig. [TBI ) indicate that E3 cannot be 
interpreted as the radio counterpart of B3/C3, since the images 
produced by El and E2 in the counter image region are located 
on the right side of B3/C3 (unlike in the case of El and E2), at 
a distance which is inconsistent with the average rms of mfinal. 
A possible scenario that can explain this result is that El and E2 
are multiple images (but we don't detect the associated counter- 
image because it is too faint) and E3 is a non-lensed radio galaxy 
that it is serendipitously lying in this region. However, as it will 
be discussed in the next section, it is also possible to explain E3 
as a lensed image produced by a source at z — 2.9 under the 
merger scenario proposed in B04. 



might also constitute a set of mirror images. To test this hypoth- 
esis, we redo the optimization using the positions of these two 
radio sources as an additional set of constraints (system 6, see 
Table |6), first assuming zcr = 2.911 and then considering zcr 
as a free parameter. A summary of the properties of the resultant 
best models (mzfix and mzfree hereafter) is presented in Table 
[10] and Table [II] 

The optimization with z free results in a predicted redshift 
for system 6 of zcr = 2.4 + 0.074, which is not consistent with 
z = 2.911. If we compare the probability of mzfix and mzfree 
to reproduce the observations (the evidence, see Table ITTb . it 
is clear that mzfree is preferred over mzfix. Note also that the 
rms difference between the observed and predicted image po- 
sitions for system 6 is three times worse for mzfix compared 
with mzfree (see Table ITOb. On the other hand, the total x 2 of 
mzfree for all the optical/IR multiply imaged systems (x 2 opti- 
cal, see Table \TT[ is just slightly worse that the corresponding 
X 2 of mfinal, which means that mfinal and mzfree are equally 
good. Note that the evidence cannot be used to compare these 
two models because the number of free parameters in each of 
them is different. 

Therefore, unless the positional uncertainty of CR1 and CR2 
could be as large as ~ 1" (which seems unlikely, since the po- 
sitional error for these sources derived from the observations is 
0.2" - 0.3"), these results indicate that the lens model favors an 
scenario in which CR1 and CR2 (if assumed to be mirror im- 
ages) are produced by a radio source that is not located at the 
same redshift as the optical arc. However, as discussed in Sect. 
13.21 CR1 and CR2 seem to be two relatively compact (~ 2") re- 
gions of an extended (~ 5") radio source. Therefore, it is possi- 
ble that this extended radio emission is a multiply imaged struc- 
ture produced by a source at z — 2.9, in which CR1 and CR2 
are not mirror images. Until the structure of this extended radio 
emission can be robustly mapped with deeper observations and 
included in the modeling process, the posibility that CR1 and 
CR2 are associated with a lensed source at redshift 2.9 remains 
open. 

As mentioned in Sect. 14.11 the multiple images of ERO B 
have a non-negligible probability of being associated with the 
radio detections El, E2 and E3. However, the probability that 
E3 and B3/C3 are associated is significantly lower than for El- 
Bl and E2-B2. For this reason, we used mfinal to check if the 
positions of El, E2 and E3 would be consistent with a set of 
three images produced by a radio source at z — 2.9. This was 
done in the following way: for each radio detection, its position 
was traced back into the source plane, and the resultant source 
was lensed again into the image plane. In this way, the model 
provided three predicted positions for each radio source, one of 
them being the position used as input. 



6. Discussion: the merger scenario 

Figure [8] shows a detail of the HST image at the position of the 
sub-mm peak, where ARC1 and the EROs are located. The con- 
tour map of the right panel indicates the total radio emission 
observed in this region, whereas the left and central panels show 
the compact and extended component. The black line correspond 
to the critical curve at z — 2.9 predicted by the lens model of the 
cluster described in Sect. [5] 

If ARC1 and the EROs are indeed multiple images produced 
by three different regions of a merger at z = 2.9, the rest of 
the merger lying between them should also be multiply-imaged. 
Hence, if the merger contains a lot of dust (as expected from a 
high-z SMG), this images would not be seen in the optical, but 
they would show up in the radio. 

Under this scenario, the fact that CR1 and CR2 do not seem 
to have optical counterparts, and lie in the region between ARC1 
and the ERO images, suggests that they could be lensed im- 
ages associated with the highly-obscured center of the merger. 
However, as discussed in Sect. 15.41 they can also be interpreted 
as a set of mirror images produced by another radio galaxy lo- 
cated at a different redshift. The real strong evidence in favor 
of the merger hypothesis comes from the extended component 
in which CR1 and CR2 seem to be embedded, which provides 
38% of the observed radio flux in the region of the sub-mm peak. 
As is clearly shown in the central panel of Fig. [8] this extended 
component constitutes a bridge of emission between the com- 
pact radio sources, as expected from the multiply imaged dust- 
obscured material of the merger located between the optical and 
NIR emitting regions. In complete agreement with this scenario, 
the shape of the total radio emission in the right panel of Fig. [8] 
follows the distribution of the optical arc and the ERO images 
remarkably well. 

In this case, we can also expect the extended emission to be 
multiply imaged like ARC1 and the EROs, producing a counter- 
image in the region between ARC1 ci and B3/C3 (see Fig.[T4li. 
Since the images in this region are de-magnified, it is likely that 
the expected radio counterpart of B3/C3, and the counter image 
associated with the extended emission, are blended. That would 
shift the peak of the total observed radio emission towards ARC1 
ci, explaining why the probability derived for E3-B3/C3 is much 
lower than what it is found for El-Bl and E2-B2. 

Finally, it is also not surprising that El-Bl and E2-B2 are not 
identified as robust (> 95%) counterpart pairs. With the cluster 
magnification generating a high resolution (distorted) view of 
the inner structure of the merger, the observed offset between El- 
Bl and E2-B2 could be real, indicating that the NIR and radio 
emissions are arising from slightly different regions inside the 
merger. 
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7. Summary and Conclusions 

SMM J04542 - 0301 is an elongated region of bright sub-mm 
emission located in the core of the clu ster MS0451.6 - 0305. It 
has been suggested dBorvs et alj|2004l) to be partially produced 
by a multiply-imaged z — 2.9 merger which contains a LBG 
(lensed as ARC1 and ARC1 ci), and two EROs (lensed as Bl, 
B2, CI, C2 and B3/C3). Given the low resolution and poor 
positional accuracy of the sub-mm map, it is very difficult to 
confirm the connection between SMM J04542 - 0301 and the 
optical/NIR lensed images. However, high resolution radio in- 
terferometric observations can help to establish this connection 
thanks to the FIR-ra dio correlation, whi c h has being found to 
be valid for SMGs dKovacs et alJ l2006t IVlahakis et ail l2007t 
llbar et al.ll2008l:lMichaIowski et al.ll2009l) . 

In a previous paper dBerciano Alba et a"L"ll2007l) . we reported 
on the detection of 1 .4 GHz radio emission coincident with this 
system using VLA archival data. Now, following a more sophis- 
ticated data reduction procedure, the previous B-array observa- 
tions have being re-reduced and combined with new high reso- 
lution A-array observations. The resultant data set has been used 
to produced two radio maps of the cluster core: 

(i) a deep (~ 10 /iJy), high resolution (~ 2") map, to compare 
with the optical/NIR images. In this map we have detected 
5 sources (El, E2, CR1, CR2 and E3) located near the 
optical/NIR lensed images, and one source (RJ) located 
in the central region of the sub-mm emission (see Fig. |4). 
Despite the fact that some of the sources have SNR < 5, the 
fact that the -Act peaks lie far from the central target region 
adds confidence to the reality of the sources. 

(ii) a 15" resolution map, to compare with SMM J04542 - 0301 
(see Fig. [Toll. 

The source RJ is the brightest and most extended of these 
six detections, with no obvious counterpart in the optical/NIR. 
Although its extended morphology could be interpreted as 
an AGN signature, the degeneracy between redshift and dust 
temperature allows the observed Ssso/jm/S' 1.4 GHz at the position 
of RJ to be consistent with a starforming SMG that follows 
the FIR-radio correlation. In any case, the low resolution of 
the 850 yum map makes impossible to determine if RJ is re- 
ally contributing to the observed sub-mm emission in this region. 

On the other hand, the other five radio detections constitute 
a strong obs ervational evidence in favor of the merger scenario 
proposed bv lBorvs et al.l d2004l) . The evidence that supports this 
conclusion can be summarized as follows: 

- if the optical arc and the EROs are multiple images produced 
by three different regions of a merger at z — 2.9, the dust 
obscured material between these regions is expected to be 
lensed in the same way, but only visible at radio and sub-mm. 
In agreement with this scenario, the shape of the radio emis- 
sion observed at the sub-mm peak follows the distribution of 
the optical arc and the ERO images remarkably well (Fig. 

— CR1 and CR2 do not have any optical/NIR counterpart. 
However, they seem to constitute two relatively compact 
emitting regions embedded in a ~ 5" extended source 
located between El and E2. The presence of this extended 
component (which contributes 38% of the total radio flux 



in the sub-mm peak) can only be explained if it is being 
produced by the dust obscured lensed material in the center 
of the merger. 

- The compact sources El and E2 have a high probability 
of being associated with the images of ERO B, but this 
probability is considerably lower for E3-B3/C3. A plausible 
explanation of this result can be provided by the merger sce- 
nario, in which the lensed extended component is expected 
to have a counter image which would be blended with the 
radio emission associated with E3/C3. This would shift the 
position of the total observed peak, lowering the probability 
of identifying E3 and B3/C3 as counterparts. 

- The fact that Bl, B2 and B3/C3 are not robust (> 95%) 
counterparts of El, E2 and E3, indicates that the radio 
and NIR emission is being produced at slightly different 
positions in the sources plane, which can be distinguished 
thanks to the enhanced resolution provided by the lensing 
magnification. 

- If the 15" resolution radio map is produced after the CC 
model of RJ is subtracted from the data, its morphology turns 
out to be remarkably similar to the sub-mm map (right panel 
of Fig. [Toll. This result strongly suggests that the sources El, 
E2, CR1, CR2 and E3 are associated with the extended sub- 
mm emission, allowing us to establish a direct link between 
SMM J04542 - 0301 and the merger. 

In this paper we have also presented a new lens model of the 
cluster MS045 1 .6 - 0305, produc ed with the latest version of the 
lenstool code (lJullo et al.l 1200 7), in which the positions of the 
ERO images and the arc knots have been used as constrains. The 
redshifts of the EROs predicted by this mo del are zb = 2.93 + 
0.13 and zc = 2.85 ±0.06. in agreement with lBorvs et al.ld2004l) . 

Including the positions of CR1 and CR2 as mirror images 
in the model significantly degrades the overall rms of the opti- 
mization. However, since these two sources are embedded in an 
extended source (which cannot be included in the lens model due 
to the current limitations of the lenstool code), it is very diffi- 
cult to establish a one-to-one identification of structures in the 
image plane. Therefore, the possibility that CR1 and CR2 are 
lensed images produced by the dust-obscured core of the merger 
remains open. 



8. Future prospects 

Radio observations of SMGs have been traditionally used to 
identify their optical/NI R counterparts, but also to provide es- 
timat es of their sizes dChapman et al.l 12004: Bigg s & Ivisonl 
120081) . On the other hand, the magnification bias due to the grav- 
itational lensing effect produced by clusters of galaxies has been 
used to increase the detection rate in sub-mm su rveys, and con- 
strain the faint flux end of the submm counts (e.g. lKnudsen et all 
120081) . The observations presented in this paper, however, illus- 
trate that the prospects for both radio interferometry and strong 
gravitational lensing in clusters of galaxies to study the inter- 
nal structure of sub-mm galaxies (SMGs) are also very promis- 
ing. With the improved sensitivity of e-MERLIN and EVLA, and 
the revolutionary view of the sub-mm sky that will be provided 
by ALMA, the enhanced resolution of multiply-imaged submm 
galaxies might permit us, among other things, to resolve for the 
first time their different star-forming regions, and help to assess 
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the connection between starburst and AGN processes that seem 
to coexist in these systems. 

SMM04542 - 0301 constitutes a very promising system to 
carry out this kind of studies, but there is still one fundamen- 
tal piece of information that is missing: redshift information 
for the NIR, radio and submm components. Current millime- 
ter facilities can provide redshift information by detecting (for 
instance) CO rotational lines, but the relatively narrow band- 
widths of their receivers (< 10 GHz compared with the 115 
GHz spacing between CO rotational lines) makes the search of 
multiple lines very time consuming, and a prior estimate of the 
redshift of the source is required. Fortunately, blind CO red- 
shift observations are now possible thanks to Z-Spec, a broad- 
band (~ 1 15 GHz) millimeter/sub-millimeter spectrograph that 
has been co mmissioned at the 1 0.4-m Caltech Submillimeter 
Observatory dBradford et alj |2009). Once all the sources that are 
contributing to the observed sub-mm emission of SMM04542 - 
0301 have been robustly identified, the increased bandwidth and 
sensitivity of the new receivers that have been installed on the 
30 m telescope (located on the Pico Veleta, Spain) and 
the Plateau de Bure Interferometer (located in the French Alps), 
will be well placed to characterize the physical properties of this 
system in detail. 

Finally note that, since these kind of systems consist of mul- 
tiple distorted images of the lensed object, a reliable reconstruc- 
tion of the source is essential for their interpretation. Therefore, 
to be able to follow this line of research, the technological de- 
velopment in the observational side has to be accompanied by a 
similar development in the lens modeling techniques in clusters 
of galaxies, to properly account for the extra complexities asso- 
ciated with extended emission and interferometric observations. 
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Fig. 2. VLA 1.4 GHz maps centered in the core of the cluster MS0451.6 - 0305. The maps were produced using the visibilities 
of the central 10 kA region of the UV-plane only, to enhance the source visibility. The facets used during the selfcalibration of the 
cluster data are indicted with boxes, while the circle represents the VLA FWHM primary beam (0.25°. at 1.4 GHz). The central 
bold box corresponds to the target region mapped in Figj4] The brightest sources in the field are indicated with labels. Top Panel: 
1° x 1° map produced with the A-array observations. Each box is 8.4' x 8.4'. Bottom Panel: 1.5° x 1.5° map produced with the 
B-array observations. The big boxes (used in the first step of the selfcalibration procedure) are 27.8' x 27.8' in size, whereas the 
small boxes (used in the second step) are 6.8' x 6.8'. 
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Fig. 4. VLA 1.4 GHz map derived from the A+B-array data set using R = 0.6. The greyscale has units of Jy beam , and the 
corresponding contours (white) are drawn at -4, 4, 5, 6, 8, 10, 14, 18, 28, 48 & 78 times the Icr noise level of 10.34 y uJy beam" 1 . 
The bl ack lines (4, 6, 7, 9, 10, 11 & 11.5 mJy beatrTl) correspond to the SCUBA 850 /mi contour map published in iBorvs et al.l 
(12004 . A comparison of the beam size of both contour maps is shown in the bottom right corner: 2.14" x 1.71" at a position angle 
of -1°.43 for the radio map, and 15" x 15" for the submm map. 
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Fig. 5. Detail of the six radio sources located within the sub-mm emission region in Fig. |4] The white radio contours are drawn 
at 2, 3, 4, 5, 6 & 8 times the Icr noise level of 10.34 //Jy beam -1 . The radio peak positions (Spk) listed in Table[T]are indicated 
with white plus-symbols. The green crosses correspond to the positions of the clean components in each source model. Orange 
plus-symbols and ellipses represent the center and FWHM of the best fitting Gaussian for each source. A scaled version of the radio 
beam (2. 14" x 1 .7 1 " at a position angle of - 1 .43) is included in the bottom right panel. 
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Fig. 7. VLA 1 .4 GHz maps obtained by tapering the A+B-array data set with a Gaussian of FWHM = 30 kA. The crosses indicate 
the positions of the radio detections identified in Fig. [4] Panel (A) includes all the emission between El and E2. In panel (B) the CC 
models of El and E2 were subtracted from the data before the tapering and cleaning procedure. Panel (C) shows what is left when 
CR1 and CR2 are also subtracted from the data. The gray-scale has units of Jy beam 1 , and the corresponding contours are drawn at 
2,3,4,5 & 6 times the lcr noise level of 14 //Jy beam -1 . The dimensions of the beam are 5.27" x 5.12" at a position angle of -88°. 25 
(panel (C), top right corner). 
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Fig. 8. Detail of the HST image in the arc region, including the critical curve at z — 2.9 1 1 predicted b y the lens model de scribed in 
Sect. [5] (black line). The red circles indicate the positions of the multiply imaged EROs mentioned in lBorvs et alj (|2004), whereas 
the blue circles correspond to the compact radio detections indentified in this region. The r adius of the circles indicate the estimated 
positional errors: ~ 0.20" for the EROs (astrometric error for the standard stars quoted in Taka ta et al. 2003). between 0.1" - 0.2" 
for the radio sources (see Table[TJ and ~ 0.08" (half a pixel) for the optical sources. The contours correspond to the ~ 2" resolution 
radio map presented in fig|4](left panel), and two of the ~ 5" resolution radio maps presented Fig. [7] (central and right panel). 
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Fig. 10. Detail of the HST image of the cluster center, with the sub-mm 850 /im contour map superimposed in black. The white 
lines correspond to the VLA 1.4 GHz contour map obtained after tapering the A+B-array data to match the resolution in sub-mm. 
Radio contours are drawn at 2, 3, 4, 5, 6 & 7 times the lcr noise level of 33.93 /iJy beam -1 . Sub-mm contours are drawn at 4, 6, 7, 
9, 10, 11 & 11.5 mJy beam -1 . The white crosses indicate the position of the radio sources identified within the sub-mm emission 
region. The beam size (15" x 15") is indicated in the top right corner. Left panel: map produced including all radio detections. 
Right panel: ma p produced after su btracting the CC model of RJ from the data. Indicated with diamonds are the positions of 3 
EROs reported in lTakata et alj d2003l) that could be contributing to the observed sub-mm emission. 
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Fig. 12. Detail of the HST image of the cluster core, indicating all the galaxies that have been included in the lens model. The 
sizes of the ellipses correspond to the morphological parameters listed in Table The lines represent the critical curves (black) and 
caustics (white) predicted by the lens model at z = 2.91 1. 





Fig. 15. Analisis of the lensing nature of El, E2 and E3 with the lens model mfinal. The orange circles indicate the image positions 
predicted after tracing each of the three radio sources independently into the source plane and lens them back into the image plane. 
The radius of the circles correspond to 0.3". The rest of the labels are the same used in Fig. [8] 



